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SUMMARY 

This analysis calculated the mass moment of inertia of a non- 
viscous fluid in a slowly rotating rectangular tank. Given the 
dimensions of the tank in the x, y, and z coordinates, the axis of 
rotation, the percentage of the tank occupied by the fluid, and 
angle of rotation, an algorithm was written that could calculate 
the mass moment of inertia of the fluid. While not included in 
this paper, the change in the mass moment of inertia of the fluid 
could then be used to calculate the force exerted by the fluid on 
the container wall. 


ANALYSIS 

Depending on the dimensions of the tank, the amount of fluid 
in the tank, and the angle of rotation, the resulting shape of the 
fluid can be broken down into simple geometries. Since the fluid 
was assumed to be incompressible, non-viscous, and the rotational 
velocity small, the shape of the fluid was represented by 
rectangular and triangular prisms. The most simple fluid shapes 
were represented by a single rectangular or a single triangular 
prism, while the most complex geometry was represented by two 
rectangular and one triangular prism (See Fig. 1 and 2). The mass 
moment of inertia ("Iz" in this coordinate system) of each prism 
was first calculated relative to its centroid. Then the parallel 
axis theorem was used to calculate the Iz of the entire volume at 
either the fluid’s center of mass or the container’s center of 
rotation . 

The mass moment of inertia of a rectangular prism with respect 

to its centroid and the axis shown was calculated by the following 

equation:^ (NASA-CR-197777) THE CALCULATION ^^^^N95-2 7165 
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where : 

= Moment of inertia with respect to the z-axis. 
a = Length of rectangular prism in x-axis. 
b = Length of rectangular prism in y-axis. 
m = Mass of prism. 

The mass moment of inertia of a triangular prism with respect to 
its centroid and the axis shown was calculated by the following 
equations : 


T O 


The parallel axis theorem:^ 




where : 

X = Distance in x-direction from centroidal to arbitrary axis, 
y = Distance in y-direction from centroidal to arbitrary axis. 


For comparison, the effective moment of inertia of fluid was 
calculated by the following formula- 

— ^-1- +2. 510 (tanh-;:^ +0.0045) ( — ^— ) 

1+Ti ^^1 1+Ti 

where : 

Ip^ = Effective moment of inertia of fuel about z-axis. 

= Moment of inertia of solidified fuel about z-axis. 

= Iz when theta equals zero, 
rl = b/a = tank aspect ratio in xy-plane. 


RESULTS 


As the partially filled rectangular tank shown in figures 1 
or 2 rotates about its origin, the Iz of the fluid will change. 
Program MomentOf Inertia (See Appendix B) calculated the Iz of the 
fluid relative to the centroid of the fluid. In tables one through 



three, three different tank dimensions (all with unit depth) are 
shown: 1 x 2, 1 x 4, and 1x8, respectively. With each tank 
dimension, the mass moment of inertia for three different fluid 
volumes were tabulated. Only zero through 90 degrees were 
calculated since the mass moment of inertia for 90 through 180 
degrees are mirror image of the shown data. The data then repeats 
every 180 degrees. This data is also shown graphically in figures 
3 through 6 . 

As expected, the Iz of tanks with aspect ratio of 1/2 did not 
change significantly relative to tank rotation at any fluid level. 
This was due to the proximity of the aspect ratio to unity. With 
smaller (or greater) aspect ratios, the change in Iz increased 
significantly. For 50% volume, there was a 112% increase in Iz/rho 
for a/b = 1/2; while for a/b = 1/8, the change was 278%. 
Decreasing fluid volume also increased the change in Iz as the tank 
rotated. For a/b = 1/8, and the fluid volume was 80 percent, the 
change in Iz was 54%. For 20 percent fluid volume and the same 
aspect ratio, the change in Iz was 1,700%. 

Program Momentl (See Appendix B) calculated the Iz of the 
fluid relative to the center of rotation, which in this case was 
the origin. The calculations were very similar to the previous, 
except that the center of rotation was used as the axis instead of 
the C.G. of the fluid. Tables 3 through 6 shown the values 
calculated for 3 different aspect ratios and 3 different fluid 
volumes for each aspect ratio. These values are shown graphically 
in figures 7 thru 9. 

As compared to the Iz relative to the C.G. of the fluid, the 
Iz relative to the center of rotation was greater , given the same 
aspect ratio and fluid volume. The greatest increase was the 20 
percent fluid volume and aspect ratio = 1/2. The smaller 
percentage filled containers had greater increases than the higher 
percentage filled containers due to the greater change in the 
distance between the C.G. of the fluid and the center of rotation. 
For example, while the Iz of zero degrees, aspect ratio = 1/5, and 
80 percent fluid volume changed 17% when the axis changed from C.G. 
of fluid to center of rotation, the Iz for the same angle and 



aspect ratio, but only 20% fluid volume, changed 1,873%. 

The change in the mass moment of inertia also increased with 
decreasing aspect ratio and decreasing fluid volume, however, the 
percent change was not as dramatic. The change in Iz for aspect 
ratio = 1/2 and 20% fluid volume relative to the C.G. was 248%, 
while the change in Iz relative to center of rotation was reduced 
to 48%. 

The Iz of a tank with decreasing fluid level was also 
calculated by altering program MomentCG . The fluid level started 
at 100% at 0 degrees and decreased to 0% after 360 degrees . The 
results for three different aspect ratios are plotted in Figures 7 
thru 9. This data demonstrated that even with an aspect ratio = 
1/2, the change in Iz can be significant when the fluid level 
decreased. With even smaller aspect ratios, such as wing fuel 
tanks, the change in Iz was even greater. The "jumps" in the value 
of Iz corresponds to the a large change in the geometry of the 
fluid . 

For comparison, the effective moment of inertia for various 
tank aspect ratios were also calculated. The effective moment of 
inertia is the moment of inertia of an equivalent mechanical 
system. However, the effective moment applies to small angular 
displacements only, which was quite different from the case that 
was analyzed here. These values are plotted in Figures 10 thru 12. 
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APPENDIX A 


TABLES AND FIGURES 
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FIGURE I : TANK CONFIGURATION 




FIGURE 2: MOST SIMPLE GEOMETRY. 



FIGURE 3: MOST COMPLEX GEOMETRY. 


Aspect Ratio 


1/2 


Izz/rho 


Theta 

20% Full 

50% Full 

80% Full 

0.0 

.0387 

.1667 

.4747 

10.0 

.0386 

.1677 

.4766 

20.0 

.0385 

.1712 

.4829 

30.0 

.0380 

.1780 

.4953 

40.0 

.0361 

.1903 

.5178 

50.0 

.0361 

.2125 

.5440 

60.0 

.0411 

.2552 

.5703 

70.0 

.0553 

.3154 

.5989 

80.0 

.1019 

.3454 

.6183 

90.0 

.1347 

.3542 

.6187 

% Change 

248 

112 

30 


Table 1: Iz relative to fluid center of gravity for tank 

dimensions: a = 1.00, b = 2.00, and c = 1.00. 


Aspect Ratio = 1/5 


Izz/rho 


Theta 

20% Full 

50% Full 

80% Full 

0.0 

.1093 

.8333 

2.9973 

10.0 

.1101 

.8358 

3 .0014 

20.0 

.1126 

.8439 

3.0147 

30.0 

.1173 

.8599 

3 .0410 

40.0 

.1256 

.8891 

3.0894 

50.0 

.1400 

.9450 

3.1825 

60.0 

.1642 

1 .0651 

3.3819 

70.0 

.2213 

1 .3867 

3.6922 

80.0 

.4158 

2.3042 

4.1547 

90.0 

1 .0693 

2.7083 

4 .4373 

% Change 

878 

225 

48 


Table 2: Iz relative to fluid center of gravity for tank 

dimensions: a = 1.00, b = 4.00, and c = 1.00. 



Aspect Ratio 


1/8 


Izz/rho 


Theta 

20% Full 

50% Full 

80% Full 

0.0 

.4747 

5.6667 

22.378 

10.0 

.4766 

5.6718 

22 .387 

20.0 

.4829 

5.6885 

22.467 

30.0 

.4953 

5.7216 

22.467 

40.0 

.5180 

5.7826 

22.566 

50.0 

.5610 

5.9000 

22.755 

60.0 

.6519 

6.1581 

23.173 

70.0 

.8833 

6 .8869 

24 .368 

80.0 

1 .6633 

10.522 

28.232 

90.0 

8.5387 

21 .417 

34 .475 

% Change 

1699 

278 

54 

Table 3: Iz 

relative 

to fluid center 

of gravity 

dimensions: a 

= 1.00, b 

= 8.00, and c = 1 . 

,00. 


Aspect Ratio = 1/2 


Izz/rho 


Theta 

20% Full 

50% Full 

80% Full 

0.0 

.2947 

.4167 

.5387 

10.0 

.2931 

.4167 

.5390 

20.0 

.2880 

.4167 

.5453 

30.0 

.2780 

.4167 

.5553 

40.0 

.2596 

.4167 

.5737 

50 .0 

.2387 

.4167 

.5946 

60.0 

.2187 

.4167 

.6147 

70.0 

.1986 

.4167 

.6347 

80.0 

.1924 

.4167 

.6409 

90.0 

.1987 

.4167 

.6347 

% Change 

34 

0 

18 


Table 4: Iz relative to center of rotation for tank dimensions: 
a = 1.00, b = 2.00, and c = 1.00. 



Aspect Ratio = 1/5 


Iz/rho 


Theta 

20% Full 

50% Full 

80% Full 

0.0 

2.1573 

2.8333 

3.5093 

10.0 

2 .1542 

2 .8333 

3.5124 

20.0 

2.1441 

2.8333 

3.5226 

30.0 

2.1240 

2.8333 

3.5427 

40.0 

2 .0869 

2.8333 

3.5797 

50.0 

2 .0153 

2.8333 

3.6514 

60.0 

1 .8067 

2.8333 

3.8060 

70 .0 

1 .6238 

2 .8333 

4 .0428 

80.0 

1 .2927 

2.8333 

4 .3740 

90.0 

1 .1973 

2 .8333 

4 .4693 

H Change 

44 

0 

27 


Table 5: Iz relative to center of rotation for tank dimensions: 
a = 1.00, b = 4.00, and c = 1.00. 


Aspect Ratio = 1/8 


Iz/rho 


Theta 

20% Full 

50% Full 

80% Full 

0.0 

16.859 

21 .667 

26.475 

10.0 

16.852 

21 .667 

26.481 

20.0 

16.832 

21 .667 

26.501 

30.0 

16.792 

21 .667 

26 .541 

40.0 

16.718 

21 .667 

26.616 

50.0 

16 .575 

21 .667 

26 .759 

60.0 

16.259 

21 .667 

27.075 

70.0 

15 .349 

21 .667 

27.984 

80.0 

12.413 

21 .667 

30.920 

90.0 

8.7947 

21 .667 

34 .539 

% Change 

48 

0 

31 


Table 6: Iz relative to center of rotation for tank dimensions: 
a = 1.00, b = 8.00, and c = 1.00. 
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Figure 4: Iz/rho Relative to Center of Rotation for 1x2 Tank 
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Figure 6: Iz/rho Relative to Center of Rotation for 1 x 6 Tank 
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APPENDIX B 
PROGRAM LISTINGS 



'ijraT. , a: 

“r ^ i :i::u:ite= :n? fTiass fiicr:?:: o* ;nert:a cr tn? f-iiic 'n a 

'e'rr^'f :c r^n^ := d^rriensioneL &= 5 

th? r-dir. ': ■ y-d'-iuD;, anc 'C tne :-C!r(into tha zizi . 

lenttr or the tanK •= alsc tr;? cer, ter rotatijn and cr'r-’'-e5Dond=- ‘l ■ 

“ne tank '“otatas lounterclcckwiee anc tn a 
This- vers ! on was ^ast itiodi tied on Z/5/9C. 


oriQtn oi tne coordinates, 
aroie is aetineo as theta. 


, ne j 
} 


PasPrinter; 


"fiD = 1,0; 


ar 


Number = ar^^ay 

L I . .90] of real ; 

ou tFi I e ;tex 

t; {file DOinter for output) 

outFileNamie rs: 

r f n g ] 1 5 j ; 

outDutl iNumDer; 


DutDiit2;Numbe?'; 


Dr in ter; char ; 


i. 

{length of tank) 

b, 

{heioht of tank) 

% 

{width cf tank) 

angle, 

(the angle of rotation in deorees) 

area, 

{total area of fluid in plane) 

I ncremen t . 

(Increment by which angle increases) 

aTr langle, 

{Max area of tnanqluar shape) 

bTr s angl e , 

(Max ares of triangle and rectangle combined 

percent , 

(Fluid level as a percent of height ’’b”) 

mass. 

(Mass of fluid) 

mcmient , 

(Mass mcment of inertia cf entire shape) 

beta. 

(2 - theta) 

tanqentBeta, 

(tandeva.'} 

tangentTheta, 

(tan( theta)) 

theta, 

(Angle cf rotation CCW in radians) 

zero;rea 1 ; 


form: integer; 

(defines the geoirietry) 






''oced^jre r:ounterRctat?(‘yan ax, ay, axl, ayl:reai); 

This procedure converts local Cord, to olobal Cord, given angle theta. 
?qin 

ax 1 = (axl 4 cosuheta). - (ay! ^ sinCtheta)); 

av := (axl * sin(theta)) ^ (avl * ccs‘- theta)); 

id; 


jnetion IBarTrl (var mini ,ta,tb:real) :rea! ; 

*his function will calculate the mass m&Ttent ot inertia of a triangular } 
jrism given its mass, base, and height. } 

?gln 

IBarTri := (mTri/13) * (scr(ta) + sqr(tb)); 

id; 


;nct i on IBarRectCvar mRect ,ra,rb:real ) :real ; 

'his function will calculate the mass moment of inertia of a rectangular ) 
-Tism given its mass, length, and height, ) 

^Qin 

IBarRect := (mRect/12) * (sqr(ra) + sqr(rb)); 
d! 


ocedure Centroid(var 


xBar,yBar , 
al ,xl ,yl , 

vv y O 

cL*. , , 


a3,x3,v3;real ) ; 


ORIGINAL PAGE IS 
OF POOR QUALITY 



th? X 


t b.s c r r 


‘.hr 


;.AreaL, 
yAreajrBtfl ; 

?Gin 

xArea j= * xl) t u2 x2) va3 ^ ac'; 

vArea ;= (al ^ vl) ^ (a2 * >2; ^ (a^ * vi'j 

x3ar ;= xArea/^area; 
vBar 1 = yArea/area: 

Id; 


jncticn WhatTvpe: integer; 

•^Th i s function deteriTiines uihat shape the fluid is in, T/pe 1 = tn angle; *) 

ftyDE 2 = triangle aro rectanoie; and type 6 = one tri angle anc two 
^rectangles. 


aArea, alpha, bArea ;real ; 

1 

eTri angle := 0.5 * sqp(a) * tangentlheta; 

aArea := aTri angle * a + (b - (a ♦ tangentlheta)/;: 

aipha »” arc. an ( b/’ a ■■ ; 

if tangentlreta O O.C then 

begin 


bTri angle ;= abs((0 

.5 ♦ sgr(b'' 

)/t srgentTheta) ; 

bArea := bTri angle 

+ abs((b * 

(a - (b/tangentTheta 

end; 



if (angle = 0.0) then 


a rectangle, } 

WhatType ;= 4 



else if (angle = 9C.0) 

then (It'S 

still a rectangle.) 

WhatType := 5 



else if theta <= aloha 

then 



) ) ) ; 


begi n 

if (aTri angle / area; then 
WhatType := i 

else if '.'aArea ) area; and 'bArea ; area; then 
uinatType ;= 2 

else UihatType := 3; 
end 


beg.i n 

ft (bTri angle ) area) then 
lOhatType := 1 

else if (bArea > area) then 

UhatType ;= 6 

else WhatT/pe := 7; 

end; 


ocedure TypeOns; 

his procedure will calculate the mass •iicn'<ent of inertia of a triangular 

risrri. 

r 

a2, b2, h, j, ex, cysreal; 
gin 

h ;= sprt((2 * area)/tangentTheta) * tangentlheta; 
j ;= sqrt((2 ^ area)/tanGentTheta) ; 
a2 := (J/3) - (a/"2); 
b2 ;= (h/3) - (b/2); 

CounterRctate(cx, cy, a2, b2); 
mofTient ;= IBarTri(iT;ass,h,j); 


UR 

Gf 


-f*cAL 

POOR QUALITY 


""vroTliT. ! 



i ?roc?flurs 

vv ■ ' ca^cuiare :ne fTiass 

:ar De C'OKir. 

oowr in::, one criano. lar 

ir 


arealri . 


areaRectl , 


h, 

iHeiqht of Rect) 

massTri , 

{Mass of trianGcOar prismi: 

massRect, 

(Mass of Rect, prism) 

mciirientTri , 

(Moment of inertia of tri 

marientRect, {Moment of inertia of Rec 

j, Qx, gy\ 


x,y, 


xCen , 


yCen, 


>i2, x2, 


w 

A iJ , / 1-“ , 


X 4 , ^ r 


x5, x5:real 

, 

eoin 


h ;= <area 

“ atr i any] e)/a; 

J := a tanoentlheta; 

arealri := 

a:r ' anole; 

areaRectl i 

:= area - aTnanqie; 

maaslri i~ 

areaTri * : * rhoj 

massRect := 

• areaRectl * c * rho; 

Qx != a./ 3 " 

* ay z ; 

gy := h + J 

i / u “ b.y z ! 

CountarRotat 

e(x2,y2,QX,gy); 

X := 0; 


y := h/2 - 

b/'2 J 

Coun ter Rotate ( x3, x3, X ,y) ; 

Centroid'ixCen ,yCen ,areaTr i ,x2,y2,areaRect 

x4 xLen 

- x2; 

y4 ;= yCen 

- v2; 

x5 ;= xCen 

- x3; 

v5 := vCen 

•“ y5; 

mcmentTri : 

= IBarTri (massTri , a, J) t 

mccientRect 

i= 1 Bar Re c t (massRe : t , a , h ) 

moment := momientiri t mc^riertRect ; 

id; 


'ocedure Type i hree j 

'his Procedure 

will calculate the mass i 

TfList be broken 

down into one tri angular 

ar 


h , J , k , iT; , 


areaTr i , 


areaRectl , 


areaRect2, 


massTri , 

(mass of triangle) 

massRectl , 

(mass of rectangle beloi,v' 

massRect2, 

(mass of rectangle next 

momentTr I , 


iTiDiTientRectl , 


mc^TientRect2, 


tri , 

{the triangular area not 

xCen, yCen, 

(centroidal coord of the 

X 1 , y 1 , 


x2, yz;, 


x3, x3, 


x4, y4, 


x5, y5, 


x6, y6, 


X 7 , y / J 


x8, yS, 


X?, yp, 



rriZtfT:?.''! !’* :n0^T r 

anc one ^eotano-'a 


■ £,Ti. 


' sjTi re ! 


to O'' t G I n . 


:o oriqin) 


■'C , zero, zero) j 


+ (maesRect ^ (sQr(xj) + 5qr(y5))) 


and 


t of inertia o* 
two rec t anou 1 ar 


a uolufTie which) 
DrieTT:5-» } 


tr ianole) 


tilled 

fluid) 


with fluid) 





ec:n 


:= ,i d; - 

= scrt((2 ^ tr i A ■ tar^otn tThe 
= j ^ *inc 9 n:Th&:i; 

K j= a “ j I 
n ‘ hi 

areaTri := 0.5 ^ h * ji 
areaRectl ;= a * mi 
areaRect2 := n ^ K: 
massTri i= arealri ^ c 
massRectl i= areaRecti 
ma=sR9ct2 := ar9aR?ct2 
x 2 := k + j /5 - ai/ 2 : 
y2 := b/2 - \2 ^ h/3); 
Cc'unterRDtate(xl „yl ,x2,y2) 
x5 := C| 

>j ;= rri/2 - b/2; 

Counter Rot ate (xc,:/6,x5,/5) 


rhc; 

c ^ rno; 
: rho; 



:= l-;/2 - 

a/'2 : 

>•4 

:= b/2 - 

h/ ; 

Counter-Rotated 

x3^y3.x 

Cen 

tro! d(xCen 

yCen ,ar 

X? 

;= xCen - 

xi I 

>7 

;= yCen - 

y: ; 

xS 

;= xCen - 

x6; 

y'd 

:= yCen - 

y<t 1 

X? 

;= xCen - 


yp 

:= yCen - 

7 > 

momentTri := 

IBarTr i 


, areaRecti ,x6,v6,areaRect2.i : , x3) j 


momientRectl 

momientRect2 


i cmassTr i ,h . j 
IBarRect'imassRec^ 


; ■^Cmaesi r i * <;. eqru.' 
1 .a.m^-^^rriassRec tl * 


;= IBarRect(iria5sRect2,n ,k ) + (massRect2 


- 5qr(yv) 
.5ar(xS) + 
i';r(x?) t 


)); 

5qr(yS) ) / ; 
Eqr(yy)); ; 


rriOH’ient ;= momentTri + mcuTieritRectl t iriC^rrientRect^ 


nd; 


^ocedure T/peFour; 
ar 

c4;real ; 
egin 

;= percent ^ b; 
fTiDfT^ent := IBarRect(rr;aes,a,c4;' ; 

TGI 


"Dcedure TypeFi ye ; 
ar 

c5;rea] ; 
eoin 

c5 ;= percent * a; 
moment := IEarRect(ma5s,c5,b) ; 
Id; 


-Dcedure TypeSix; 

'his Procedure will calculate the mass moment of* inertia of a volume which) 
ran be broken down into one triangular and one rectangular prisrrj, } 

ir 


areaZ, 
arealri , 
areaRecti , 

(Area of 

Rect) 


h. 

{Height 

of Rect) 


massTri , 

(Hass of 

triangular 

prism) 

massRect. 

{Mass of 

Rect. pri 

sm) 

momentTri , 

(Moment 

of inertia 

of triangular prisrrj 

momentRect, 

(Moment 

of inertia 

of Rect. prism rel. 

J, ox, gx, 
tri . 

(area of 

tank' not 

filled with fluid) 


rel. to origin) 
tc origin) 
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yj.. 

y 4 :reai ; 



arealri := 0.5 ^ sqr(b) * tangentB^ta; 
areaRectl := area - arealri; 
h ;= (area - areaTri)/b; 
j ;= b * tangentBeta; 
maasTri areaTri * : ♦ rhc; 
massRect != areaRectl * c ^ ^fici 
Qx 1“ h + ~ a./ 

Qy := fc./s - b/2; 

CcunterRotate(x2,y2,9x ,gy) ; 

X := (h/2) - (a/2); 


nd; 


Cc'jnterRotateCxl ,yl ; 

Cen tr j i d( X Den »yCen, AreaTri 2, >2, areaRectl ,xl ,yl ^zerQ, 2 ero,zero) ; 
x2 1= xCen * x2j 

v2 1 = yCen - >2; 

yA := xCen - xl; 

y4 ;= yCen - yl ; 

mofTientlrl ;= IBarTr : (massTri ,b,j) (iTiassIri ^ (sqr(x3) t 5qr(y3))); 
mc«mentRect := iBarRect(mas5Re:t ,b,h) + (rriassRect ^ (sqr(x4) ^ e.qr(v4;))i 

fTiC4Tient ;= mc^rientlri + mofnentRect; 


rocedure TypeSeven; 

This Procedure will calculate tne mass moment of inertia of a volume which) 

must be broken oovjn into one triangular and two rectangular prisrrfs. 

ar 

h ^ J » k I fTi n 

areaTri , 
areaRectl , 
areaRect2, 

massTri, {mass of triangle} 

massRecti, (mass of rectangle beloc'-j triangle) 

rTiaeER9ct2, (mass of rectangle next to triangle) 

momen tTr i , 

mcmentRec tl , 

rr(C'mentRect2, 

tri, (the triangular area not filled with fluid) 

xCen. yCen, 
axl, ayl, 
ax2, ay2. 
ax 3, ay3, 

xl, x2, x5, x4, x5, xi, 

yi, y2, y3, y4, y5, yo, 

moimentl , moment2:real ; 
egin 

tri != (a b) - area; 

j ;= sqrt((2 * tr i )/tangentBeta) i 
h := j * tanqentBeta; 

k ;= fa “ j; ^ 

m }= a - h; 

areaTri ;= 0.5 * h * j; 
areaRectl := b m; 
areaRect2 := h ^ k; 
massTri := areaTri ^ c rho; 
massRectl := areaRectl c rho; 

massRectZ ;= areaRect2 ^ c * rho; 

x2 := a./2 - (2 ^ h/3); 
v2 := b/2 - (2 j/3); 

CounterRotateCxl ,yl ,x2,y2) ; 
xj := m.-'2 - a/’2; 
v5 j= 0? 
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CounterRo'cTt ^y4): 

C?n:raid(xC-er ,/CenareaTr i ,xl ,y: ^areaRecti ,x^,y^,ar?aRect2,x3,>3) ; 

axl ;= xCer - xl; 

a/I := /Cef^ - yi j 

ax2 := xCer - >:6i 

a/2 := /Can - y 6 \ 

axS := x3en - x3; 

av‘3 ;= yCen - :/3; 

momentTri ;= IBarTri (massTri ,n, j)t(fnas5-Tri (sqr(axl) ^ =qr(a/l))); 

riomentRectl := IBarRectCmassRectl ,D.rr^)T(iTi£S5RectI * (sQr(ax2)?5qr(ay2 ))> i 

mciTiBntRec t2 ;= IBarRect(ma5afiect2 ,h ,r;) + (rr»as5Rect2 * (5qr(ax3)+SQrvay3))) ; 

maTient := iriC^Tjantlr i + mofTier tRectI + rr}OfrientRect2; 


agir (rSTART ^i4;N PRDGRAM^^') 


Cl earSc'een ; 
zero := D: 




wr 1 1 9 1 n '.i " Enter d i me n s i on s 
read/ (a/,:; ; 

for 

rectanqul ar 

tank ; 

write In r'Erter water leue 
readln (percent) ; 

1 i'r- 

a oercentage 

of ’ 

while ((percent ;= 1) or ( 
beoin 

perce 

nt \= 0)) do 


wri te]n(/;.0 < Water Leuel 

read] n(pe''cen t) ; 

< :,0. Try 

again 


end; 

wri tein(’''Ente:' angle theta increrrient. ' ); 
readln( increirient) ; 

wri teInC'Enter output tile ra/ne. '); 

read] n(ou tFi ] eNarrie) ; 

angle ;= O.’J; 

area ;= s. * t ^ percent; 

masa := a * Eercent * b * 

rewr i te(:'Ot^ - e ^ootFi 1 eName) ; 

ovr ^ t?ln(outFi " ,a;b;4) ; 

wri te]n(oi::Fi ie/B = ■',b:6:4;; 

wri telnCc'utFi ! e/C = Sc:8;4); 

tvri tel n( cu tFi ] e/' Percent full = %percent:B:4) ; 

v/r i te } n(outFi : e /Type Theta Moment") ; 

v\»h i I e ^ an Q I e <= ?0 . 0 ) qc- 

Degin 

(TronTient != 0,0; 

beta ;= (FO - angle)/360)*(2 * pi); 

theta ;= ■ angl e/5'60)^(2 * pi); 

tangent Bet a ;= sin(beta)/cD5(beta) ; 
tanoentTneta ;= sln(theta)/CD5( theta) ; 
form ;= WhatT/pe; 
case ( forr.) of 
1 iTypeOne; 

2;T/peTwo; 

3:T/peThree ; 

4;T/peFour ; 

5;TypeFiue; 

i:T/peSix; 

7;TypeSeven; 

end; 

wr i tel n(outFi 1 e ,f orm, " /an9le:6:2!i " ", moment ;S:4) ; 

angle ;= angle + increment; 
end; 

c!ose(cutFi le'M 
0. (*END MAIN PROGRAiM#) 





^rooraifi Momentl: 

■ This program calculates the mass moment of inertia oi the fluid in a tank ) 

^relative the to the center of the tank. The tank is dimensioned as 'a' in ) 

ithe x-dir, 'b' in the y-dir(up), and 'c' in the z-dir(into the page). The ) 

icenter of the tank is also the center of rotation and corresponds to the ) 

^origin of the coordinates. The tank rotates counterclockwise and this ) 

'angle is defined as theta. This version was lase modified on 2/4/90. ) 

Jses 

PasPrinter; 

:^onst 

pi = 3.141593; 
rho = 1.0; 
ype 

Number = array El. .90] of real; 

'ar 

outFile :text; (file pointer for output) 
outFileName :stringE15]; 
outputl :Number; 
output2:Nufflber; 
printer:char ; 

b, 

c, 

area, 
angle, 
increment, 
aTri angle, 
bTri angle, 
percent, 
mass, 
moment, 
beta, 

tangentBeta, 
tangentTheta, 
theta:real ; 
form: integer; 


(length of tank) 

(height of tank) 

(width of tank) 

(the angle of rotation in degrees) 

(Increment by which angle increases) 

(Max area of triangluar shape) 

(Max area of triangle and rectangle combined) 
(Fluid level as a percent of height 'b') 

(Mass of fluid) 

(Mass moment of inertia of entire shape) 

(2 - theta) 

(tan(beta)} 

(tan(theta)) 

(Angle of rotation CCW in radians) 

(defines the geometry) 


»»«»»»»»»»«*»*»**»«»»»»»»»***»»»*»**»»»«»»*»*»**»*»*»»*»♦»*»***********•***) 


rocedure CounterRotateivar ax, ay, axl, aylrreal); 

This procedure converts local Cord, to global Cord, given angle theta.) 
egin 

ax := (axl * cos(theta)) - (ayl * sin(theta)); 

ay ;= (axl * sin(theta)) + (ayl » cos(theta)); 

nd; 


unction IBarTri(var mTri ,ta,tb:real)!real ; 

This function will calculate the mass moment of inertia of a triangular 

prism given its mass, base, and height. 

egin 

IBarTri := (mTri/IB) * (sqr(ta) + sqr(tb)); 
id; 


jnction IBarRect(var mRect,ra,rb:real)!real ; 

This function will calculate the mass moment of inertia of a rectangular 
prism given its mass, length, and height. 

?gin 

IBarRect := (mRect/12) * (sqr(ra) ♦ sqr(rb)); 
id; 


inction UhatType; integer; ... . i • , 

♦This function determines what shape the fluid is in. Type 1 = triangle; *) 

♦type 2 = triangle and rectangle; and type 3 = one triangle and two *) 
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<ar 

aArea, alpha, bArea treat; 

■eg in 

aTrianglt := 0.5 » sqr(a) » tangentTheta; 

aArea •= aTriangle ♦ (a * (b ~ (a * tangentTheta))); 

alpha arctan(b/a); 

if tangeatTheta <> 0.0 then 

begin 

bTriangle ;= ab5((0.5 » sqr(b))/tangentTheta) | 

bArea s= bTriangle + abs((b » (a - (b/tangentTheta)))); 

end; 

if tangle * 0.0) then (It's a rectangle.) 

UhatType := 4 

else if (angle = 90.0) then (It's still a rectangle.) 
UhatT/pe ;= 5 

else if theta <= alpha then 
begin 

if (aTriangle ) area) then 
WhatType := 1 

else if (aArea > area) and (bArea ) area) then 
WhatType := 2 
else UhatType ;= 3; 
end 
else 
begin 

if (bTriangle > area) then 
WhatType ;= 1 

else if (bArea > area) then 

WhatType := 6 

else WhatType ;= 7; 

end; 


rocedure TypeOne; 

This procedure will calculate the mass 

prism. 

ar 

a2, b2, h, J, cx, cyireal ; 
egin 

h := sqrt((2 » area)/tangentTheta) 
j ;= sqrt((2 » area)/tangentTheta) ; 
a2 := (J/3) - (a/2); 
b2 ;= (h./3) - (b/2); 

CounterRotate(cX) cy, a2, b2); 
moment := IBarTri(mass,h,j) + (mass 
nd; 


mtxnent of inertia of a triangular ) 

) 


tangentTheta; 


»(sqr(cx) + sqr(cy))); 


rocedure TypeTwo; 

This Procedure will calculate the mass moment of inertia of a volume 
can be broken down into one triangular and one rectangular prism. 


ar 

area2, (Area of Rect) 
h, (Height of Rect) 

massTri, (Mass of triangular prism) 
massRect, (Mass of Rect. prism) 

momentTri, (Moment of inertia of triangular prism rel . to origin) 
momentRect, (Moment of inertia of Rect. prism rel. to origin) 

Ji 9/» xbar, ybar, 
x, y, x3, y3;real; 
egin 

h := (area - atriangle)/a; 

j := a i tangentTheta; 

massTri := atri angle * c » rho; 

massRect ;= (area - atriangle) * c * rho; 

gx ;= a/3 - a/2; 

ov != h + i/3 - b/2: 


which) 

) 
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CounterRotateCxbarjybarjQXiQy) ; 

momentTri := IBarTrKmassTri ,a,j) + (massTri * (sqr(xbar) ♦ 5qp(ybar))); 

X := 0; 

y ;= h/2 - b/2; 

CounterRotate(x3,y3,x,y); 

momentRect := lBarRect(mas5Rect,a,h) + (massRect » <sqr(x3) ♦ sqr(y3))); 
moment ;= momentTri + momentRect; 

?nd; 


Procedure TypeThree; u i.\ 

iThis Procedure will calculate the mass moment of inertia of a volume whichf 

imust be broken down into one triangular and two rectangular prisms. ) 

»ar 

h| Jf kj m^ 

massTri, (mass of triangle) 

massRecU, (mass of rectangle below triangle) 

mas5Rect2, (mass of rectangle next to triangle) 

momentTri , 

momentRectl , 

momentRect2, 

tri, (the triangular area not filled with fluid) 

xl, x2, x3, x4, x5, x6, 

yl, y2, y3, y4, y5, y6, 

momentl, moment2:real ; 

}egin 

tri := (a » b) - area; 
j ;= sqrt<(2 * tri)/tangentTheta) ; 

h := J » tangentTheta; 

k := a - j; 

m != b - h; 

massTri !=D.5*h»j*c» rho; 
massRectl := a » m » c » rho; 

massRect2 ;= h * k * c » rho; 

x2 := k ♦ j/3 - a/2; 
y2 := b/2 - <2 » h/3>; 

CounterRotatetxl ,yl ,x2,y2) ; 

momentTri := IBarTri(massTri ,h,j) + (massTri * (sqr(xl) ♦ sqr(yl))); 

x5 := 0; 

y5 := m/2 - h/2; 

CounterRotate(x6,y6,x5,y5); 

momentRectl := I BarRect (massRectl ,a,m)+ (massRectl * (sqr(x<4) ♦ sqr(y6)>); 

x4 ;= k/2 - a/2; 
y4 ;= b/2 - h/2; 

CounterRotate(x3,y3,x4,y4); . ...vx 

momentRect2 := IBarRect(ma5sRect2,h,k)+(massRect2 * (sqr(x3) + sqr(y3))); 

moment := momentTri ♦ momentRectl ♦ momentRect2; 

>nd; 


■rocedure TypeFour; 

■ar 

c4:real; 

legin 

c4 ;= percent » b; 

moment ;= IBarRect(mass,a,c4) + (mass » sqr((b/2) * (percent - 1))>; 

nd; 


'rocedure TypeFiwe; 
ar 

c5:real ; 

'egin 

c5 != percent » a; 

moment := IBarRect(mass,c5,b) + (mass * (sqr((a/2) * (percent 

nd; 


1 »)); 
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Procedure TypeSix; 

(This Procedure will calculate the mass moment ot inertia ot a volume which) 
(can be broken down into one triangular and one rectangular prism. ) 

uar 

area2, (Area of Rect) 

areaTri angle, 

h, (Height of Rect) 

massTri, (Mass of triangular prism) 

massRect, (Mass of Rect. prism) 

momentTri, (Moment of inertia of triangular prism rel. to origin) 

momentRect, (Moment of inertia of Rect. prism rel. to origin) 

J> 9x> 9y» xbar, ybar, 

X, y, x3, y3:real; 
begin 

areaTriangle ;= 0.5 » sqr(b) * tangentBetaj 

h := (area - areaTriangle)/b; 

j := b « tangentBeta; 

massTri := areaTriangle * c * rho; 

massRect ;= (area - areaTriangle) » c » rho; 

gx := h + j/3 - a/2; 

gy ;= b/3 - b/2; 

CounterRotate(xbar,ybar,gx,gy) ; 

momentTri := IBarTri(massTri ,b,j) + (massTri * (sqr(xbar) + sqr(ybar))>; 

X := (h/2) - (a/2); 
y := 0.0; 

CounterRotate(x3,y3,x ,y> ; 

momentRect := iBarRect(massRect,b,h) + (massRect * (sqr(x3) + sqr(y3))); 
moment := momentTri * momentRect; 

end; 


Procedure TypeSeven; 

(This Procedure will calculate the mass moment of inertia of a volume which) 

(must be broken down into one triangular and two rectangular prisms. ) 

var 

h , J , k , m, 

massTri, (mass of triangle) 

massRecU, (mass of rectangle below triangle) 

massRect2, (mass of rectangle next to triangle) 

momentTri , 

momentRecU , 

fflomentRect2, 

tri, (the triangular area not filled with fluid) 

xl, x2, x3, x4, x5, x6, 

yl, y2, y3, y4, y5, y6, 

momentl, moment2:real ; 
oegin 

tri ;= (a * b) - area; 
j := sqrt((2 * tri)/tangentBeta); 

h := J » tangentBeta; 

k := b - j; 

m := a - h; 

massTri := 0.5 » h » j * c * rho; 

massRectl := b » m » c * rho; 

massRect2 := h » k » c » rho; 

x2 := a/2 - (2 * h/3); 

y2 ;= b/2 - (2 » j/3); 

CounterRotate(xl,yl,x2,y2); DftirilKiAi n 

momentTri := IBarTri(massTri ,h,j)+(massTri * (sqr(xl) ♦ sqr(yl))); __ "AGE IS 

x5 ;= m/2 - a/2; 0^ POOR QUAUTY 

y5 ;= 0; 

CounterRotate(x(i,yd,x5,y5) ; 

momentRect! := IBarRect(massRectl ,b,m)+(massRectl » (sqr(x^) ♦ sqr(y^))); 

x4 ;= a/2 - h/2; 
y4 := k/2 - h/2; 

CounterRotate(x3,y3,x4,y4) ; 

momentRect2 := IBarRect(ffiassRect2,h,k)+(roassRect2 * (sqr(x3) ♦ sqr(y3))); 

moment ;= momentTri * momentRect! + momentRect2; 



end; 


{«»«>»»««»»»»»«»»»»»»«»«»»«»»»»>»»»»«»»««»»»«»»»»»»»»»«»»»»»»»»»»««»»»»«»>«} 
Begin (♦START WIN PROGRAM#) 

Cl ear Screen; 

writein ('Enter dimensions for rectangular tank:'); 

readin (a,b,c); 

writein ('Enter water level as a percentage of ’b*.'); 
readin (percent); 

while ((percent >= 1) or (percent (= 0)) do 
begin 

writeln('O.D ( Water Level < 1.0. Try again.'); 
readln(percent); 

end; 

wri te1n('Enter angle theta increment.'); 

readln(increment); 

writeln('Enter output file name.'); 

readln(outFileName); 

angle := 0.0; 

area := a ♦ b ♦ percent; 

mass := a ♦ percent #b#c; 

rewr i te(outFi 1e,outFi 1 eName) ; 

writeln(outFi1e,'A = ',a:8;4); 

writeln(outFile,'B = ',b:8:4); 

wri teln(outFi le,'C = ',c:8:4); 

writeln(outFile, 'Percent full = ',percent:8:4); 

wri teln(outFi le, 'Type Theta Moment'); 

while (angle (= PO.O) do 

begin 

moment := 0.0; 

beta := ((?0 - angle)/360)#(2 ♦ pi); 

theta := (ang1e/3^0)#(2 ♦pi); 

tangentBeta := sin(beta)/cos(beta); 
tangentTheta ;= si n(theta)/cos( theta); 
form := UhatType; 
case (form) of 
1 :TypeOne; 

2:TypeTwo; 

3:TypeThree; 

4:TypeFour; 

5:TypeFive; 

6:TypeSix; 

7:TypeSeven; 

end; 

writeln(outFi1e,form, ' ',angle:6:2, ' ', moment :8;4); 

angle := angle ♦ increment; 

end; 

c1ose(outFi le); 

Ind. (♦END WIN PRDGRf^#) 



APPENDIX B 


PROGRAM LISTINGS 




('^urTiu'c' ~ 3.r“5.y 


Dut- : . &Nar*ir :e 

Ow^Dwtl jNlirD®^ f 
jutLL:2:Num:e " i 


i e do: n:t" tor outcj:.' 


^* ] p r t h C-'* t s. p. K } 

■ hticnt ot tarf,) 

{wiath ct tank) 

(tht anoiE 0 - ^r-Tation :r c! 8 Q^-&ea> 


'iC! e : r,creasp = : 


.:c:a; area -'jtc ir u 
wr.cr anoU 
{-ax a^?a o- t^;aro.-ar = 
{"'a:: :* :^;ar:'t sr: 

■“'u^c as a Denat: 


{ r r '’rC-T- a r* 0* I* 


tarq?r:::e:a, aan^ra:a.»i 

:aror' ; tan(tn?:a.O 


^ 2 

V M r V d f 


\:^nciir Ot !"o:atJcn :r *ac:;ana; 


-orTij irteoerj .oatinas th? qeonjetry.^ 


'oosdu’s Coun: 5 “Rc:at?(^<ar ax, ay, axl, a/ijraai^j 
Fnis otoredure converrs =o:a: Corb. to OiO^al Cord, qi 


ile thetc 


ax ;= ;axl t C05(theta)) - Cayl t sinaheta;^i 

av := (axl t- sin(theta)) + (ay! * cos’: :heta);: 


unction I 3 arTri(yar mTri ,ta,tb:real):rea: ; 

^nis tunotion uHI cakuiate the masa moTent o- Inertia o-f a trianoular 

on^m Q.ven its mass, base, and height, 

eoln 

IBarTri := (mTri/ 19 ) * (spr(ta) + sqr(tb))r 

nd: 


jnc t ! on I BarRec t ( y ar mRec t , ra , ro : r e al ) ; r eal ; 

"hiS function wi]] calculate the mass mccrient of inertia cf a rectanoular 

orism; giuen its mass, length, and height, 

egin 

IBarRect ;= ^mRect/lZ) '* (sqr(ra) ^ sqr(rb;); 


' oc 3 du r e Ce n t r o i d ( u ar ;< Bar , y Bar , 

al,xl,yl. 
5.2,x2,y2, 
a3,x3,y3:real ) ; 
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This DrC’Crdur? 
:? -uiC. 
ar 


ra'CuUte the x an: :::rd tne cervc^d or tne 


yf-'r E-a « " f c ‘ I 


;eci n 


;= ^al ^ xl' 
yi^rea *= ’ai vi) 
xBa** := xAr^pa/areai 
yBa" 1 = y^Hrea/a^ea: 


.O.Z ^ xd.' - 
xd 2 * > 2 ) U 3 


?nd; 


•uncticn lyhitlyoeijntcQer; 

♦Thia runction determines ^nat anaDe tne r'uid 


. >oe 


t‘ ; ano. 


^tyD& 2 = triangle an: reotangie: and t^ne 3 = one t" angle and two 

*rec tan dies. 


iRr e a , a i d n a , oRr e a ; e a m 


j e 0 ! n 


a ■ f*' : apQ ! e 
a- "" e c I “ 
a ^ : L. a ; “ 


:= 0,5 * acr(a) * tanoentiheta; 


a'";ang.e * a ^ 

1 i “ ar 1 1 an 0/' a, • \ 
anoent^nete I- T.G ther 


tanoent ineta)!^' ; 


: an: 


* Urg; r = ^ ^ tnen 

"^^natTyre i- 4 


b^’- STo'e + ada^'t ^a - o/ tangert'^ne ta 
(It'a a rectangle.) 
then ilt'a a:’’’ a 'e:tano:5,. 


‘ a “ j = 


a i D n a 


wnat.yoe := ; 

S'Se ■■aRi'ea ■ area/ and iRrea 

wna:“yoe := 2 

else whatType := 3j 


■?aJ then 


enc 
becii n 


IT (blrlangie / area) then 
Uhatlype 1 

ei=® it (oRrea ) area) tnen 

WhatType := 6 

else WhatlyDe := 7; 


end; 


Vccadure TypeGne ; 

This procedure vvill calculate the mass Ticff<ent of inertia ot a triangular 
-pr 1 s<Ti. 
ar 

a2, d 2, h, j, cx, cy;rea!; 


ecin 


n ;= 5-crtui + area.'/ tangent • neta) * targentmeta; 
j ;= :-prt((2 * area)/tancentTheta; ; 
a *! i — \ J / w ) ~ c 5 l^ d ; 
b2 ;= (h/5) - ^:/2); 

CounterRctate(:;: , cy, a2, b2) i 
fTioment ;= IBarTri(rr.aeE,h,j>; 


ginal page fs 

OF POOR QUALITY 


r.d; 


e ; vpetuic; 


rocedur 



F-'Ocedure wiM 
:ar De aown 


la’cuiate the fr.a==. 
r.:c one trianaL :ar 


tj *: * '•* I n e ^ t ? 

ar,c cne reztsno’j'' a" 


3, 

r»r ‘ sjtj. 


lAih : :r ! 


area:^! , 

S'* c iPf , 

t- : 

naiERec: . 

rr;r ’''C^- t”^r ■( 


-.-siQn: c* -:ect; 

!-l£.E5 of Df 

Rect. Dr ' EfT;} 


incfneiit ot 


[re C't t r i aT; uii : ar Dr ST. 


iT!Ofr:rrf:Rect,{Mcrfri8nt 
J, QX, gv. 

icCer, 

yCen, 


ot inertia ot Rect, or is?:: re I 


re;, to o-rgin;: 
:o origin) 


x 5 , yoireah 
ruin 

n != -area - atr? ang-E ;/a; 
j ;= a tanoentThetai 
areaTr: := af angle; 

areaRec:: i- area - aTriancle; 

KiisiTn ;= area'^r’ ^ ^ 

rriasa^ert i= o.^Er^E:*:! ^ z ^ !"hD; 




uOL^'*“''Kota:t ■ 



\ .1 ^ ^ . a ' c a ■* . r w . ^ I 




yn 


-Can 

xCen 



Tia E s I r 


ac 


v4 ) .' ‘ ; 




nen: 


- Tion^ert-ec 


"ocedure TvoeTh 

ree ; 


^his Procedure 

v^i ! 1 

calculate the friasa 

Tiust 08 broKen 

do-wn 

into one triangular 

ir 

h, j, K, rr;, 

areaTr i . 
areaRectl , 
areaRe:t2, 
masE-Tri , 


of triangle) 

fTiassRectl , 

{r:as5 

of rectangle belo^v 

massReetZ, 

(lT:aES 

ot rectangle next 

mc^entTr i , 
iTicaentRectl , 
moiTientRect2, 
tri , 

(the 

triangular area no: 

xCen, yCen. 

teen t: 

'oidal coord of the 


fTiCrTTient OT inertia of 
and two re:tanc-jlar 


triangle) 
to triangle) 


tilled with ^luid) 
fluid) 


V 

i*'3, 


x4, 

v4. 



yo, 


••• 1^. 

y6, 



y*' j 


x8, 

Vo, 


>;?. 



mc^!entl , 

fTioment 


eal ; 


a 1 ‘oluiTie which} 

r\ z i c.fT, c , 
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eojn 


j= u * D) - area; 

•s s 2 r ^ ( 2 - t r i ) t an D6 r tTh eta* * 
;= .i ^ varce* :The:a; 


- — ^ I 


E"-eai^t := w’.5 t r j; 

a r a a a c t * i — a ^ rri ; 
areaRe:t2 := r ♦ k: 
fnaae':'! ;= areaTri * c rhc; 
naE5Ke:tl := areaRectl c * rnc; 

.TiaesRec:2 ;= ar9aP.ect2 c ^ rho; 

x2 := k t j/3 - a./2: 
y2 ;= b/2 - (2 * r/3); 

2 r* _ ^ t e ^ R C- ^ ^ t e ( X ^ ^ 2 « ’> 2 '•' ! 

>t j « “ ? 

p 1 1 f' * i r' ^ ^ 2 * Q -.V yb. w ^ V ^ 'f * 

x^ 2 ~ a * *i J 

v4 := d/ 2 - h/2; 
CD2rte!''KDtate(x3,>3,x4,y4'>: 

Centres d(xCe“ ,/Cen , are air i ,x! ,yl , areaRectl 


X L , yi , are aRec 1 2 r x3 , y3) ; 




: ;= Tasair- ,r .. -i-aEs: 

tRe: tl := IBarRecttinaesFe: *1 ,a.?T - 

iTiC^ientRecrC ;= I5a^RectUia=aRec:2,* ■'*. 

rfiorjent n-icne't^ri rr.Dr:er:tRect: r r 


c =bn.t'_ . 


” t £ 2 r = ;; y ■; 

* •■ ecr ^ x 9 ) 



mor-ert ;= IE:arRecttn*:aea,a.c41’ ; 
nd; 


roc e dure TypeFive: 
ar 

c5;rea] ; 
eoin 

c5 := percent ^ a; 
inoment ;= I&arRect(mas5,c5,b) ; 
id; 


''ocedure TypeSix; 
This Procedure wil] 
can be broken down 
ar 


calculate the mass 
into one triangular 


moment oi inertia oi a volume 
and one rectangular prismi. 


which} 

} 


area2. 

(Area of 

Sect) 


area'i'r i f 
areaRec tl , 
h. 

•Height 

ot Rect) 


massTr i , 

{tiiSS of 

tri angular 

or 

massRect, 

(Mass of 

Rect. pri 

sm) 

momentTri , 

(Morrient 

of Inertia 

of 

meaner tRect. 

(Moaient 

of inertia 

of 

j, gx, gy. 
tri , 

(area of 

tank net 

fill 


i STi) 

triangular prisrri rel. to origin) 
Rect. prism rel. to origin) 

ed with fluid) 





:anient5?:a; 


arear\e:t- 


. - o.'ce *“ r e c : . 


n j= \arsr - drpa;^i>/Dj 
, := D :ar:pe?:: 5 etaj 

fTiaaalri ;= araaTr^ » rhcj 

iTiassRe:: := areaRectl + c ^ ^noi 

"tV *= r-. - • - Ji/'y* 

w/^ t w* a/ ^ « 

QV J= D-“S - b/2; 

^ - a n K r. • s t a (■ V ■“■ •,. “ w -. V ■ 1 » 

• “ I - 1 

»•• w w 1 ( '. 8 . r>. w a 1 , 8 ■'. .^ i t f ' < 1 I 

, vl;eri,Arei”r ! .i.Ir-I.^reaRecti ,xl ,vl ,z?rc 42 ero,zar 


X- := ;a.er: - xl; 

:= vCar - yl ; 

rriO?T;art“' ‘ j= IBarT" ^fTa=='^r » ,b, j > ^ ^.Aasalri ^ ^aar^x3; 

fT/jTtP" ■ “■'8 “ s := ■■ Bs.''R9c * \M dz aRa I " 1 n I h ) ^ (fTiaEr-^^act ^ '.ec’* 


' a ■.-■ - □ - a i ; 0 « 


Aiu ata :?r ”a=E Ticeant Ea if^art-a a yo' 

t: t’^iEncuia^ arE two :"e:tanou i a^ dms::; 


rEioriient ; r » , 
icemen tRec'l , 
morr-entReotB, 

xCerj. yCeru 
axl, ayi, 
a X 2 4 5 » 

aXw4 dL/Oi 


.ra== OT recEanc’E ba!Dw triancie^ 
.HidEa eA rectanQir oaxt tc- t:".an^!e} 


[the tnanquiar area net tnied witn fluid) 


f 4 A w 4 A U' 4 


yl, v2, y2, y4, yo, yd, 
fTiDfT^entl , rriKnent2:r8a] \ 


tri := (a * b) ~ area; 
j := =qrt((2 ^ tri )/tar. 2 entBeta) : 

t: ;= j '* tanoentBeta; 

k ;= b - i; ‘ 

m ;= a - h; 

areal r-i ;= C,5 ♦ h ♦ j; 

araaRectl ]= b ^ 

areaStEtZ ;= h ^ k; 

maisTr? := arsaTrl * c * rho; 

rr.aoE.Reotl := areaRectl * c * rho; 

'TiaEsRectZ ;= areaRect2 * c ^ rno; 

xZ ;= el/ 2 - (2 ^ h/3): 

y2 := b/2 - '2 j/3-; 

CounterRotc* e(xi ,yl ,;;2,>2; : 
x5 := ni/2 - a/2; 
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;,ir : ; n ■ E" ang]f tri5:a !n:r5f^?n:. 

readlr( ;n:reiTi&rt; ; 

v;r; te :r- lutDLj: -:*5 nans.''-^; 

c a^ i r '•, Cta t “ ! ' ' i 


i’fE ;= a ♦ percen:: 



vr- :?]n(L^utFi ]e,"C = ",ci6;4); 


wr; te]n(Du:^i lE/^ercent tul: = ",percEnt:9 
i-vP' teln(outFi le, Typ« Theta Moment ') ; 
while (ano'e <= ?0.G) a:- 
Deg in 

mctT^ent t= 0.0; 

j-ta ;= <(?0 - anQle)/360)*(2 * Di); 

theta ;= (anol e/36D>^(2 # pi); 

tanoentBeta ;= sin(beta)/ccs(beta) ; 
lanoentlheta ;= sin( the ta)/co5( theta) ; 
Torm := Whatlype; 
case (form) of 
1 iI'/DeOne; 

2:T>-pe7wo; 
oiTypeThree ; 

4;TypeFour ; 

5;TypeFi ve ; 
d jTypeS; 

"iTvDeSeyen; 

end; 

wr f tel n(ou:’^ 1 e Jormt ' Sangl e;6;2, 
angle ;= angle + increment; 


clo5e(outFi le); 

1C. (*END ^V:^IN PRODRAtH^) 


' ,rrn3rient:S;4); 
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^rooram Momentl; 

;:Thi5 program calculates the mass moment of inertia of the fluid in a tank ) 

^relative the to the center of the tank. The tank is dimensioned as 'a' in } 

ithe x-dir, ‘b' in the y-dir(up), and ■'c' in the z-dir(into the page). The ) 

^center of the tank is also the center of rotation and corresponds to the ) 

iorigin of the coordinates. The tank rotates counterclockwise and this ) 

iangfe is defined as theta. This version was lase modified on 2/4/90. ) 

Jses 

PasPrinter; 

!^nst 

pi = 3.141593; 
rho = 1.0; 

Type 

Number = array II. .903 of real; 
tar 

outFile :text; {file pointer for output) 
outFileName :string[15]; 
outputl :Number ; 
output2:Number ; 
pr inter :char; 

A « 

b. 

c, 

area, 
angle, 
increment, 
aTri angle, 
bTri angle, 
percent, 
mass, 
moTTient, 
beta, 

tangentBeta, 
tangentTheta, 
thetarreal ; 
form: integer; 


{length of tank) 

(height of tank) 

(width of tank) 

(the angle of rotation in degrees) 

(Increment by which angle increases) 

(Max area of triangluar shape) 

(Max area of triangle and rectangle combined) 
(Fluid level as a percent of height 'b") 

(Mass of fluid) 

(Mass moment of inertia of entire shape) 

(2 “ theta) 

(tan(beta)) 

(tan(theta)) 

(Angle of rotation CCW in radians) 

(defines the geometry) 




'rocedure CounterRotateCvar ax, ay, axl, aylrreal); 

wThis procedure converts local Cord, to global Cord, given angle theta.) 
:^egin 

ax ;= (axl * cos(theta)) - (ayl « sin(theta)); 
ay := (axl * sin(theta)) + (ayl » cos(theta))j 

?nd; 


‘unction IBarTr i(var mTri ,ta,tb:real):real ; 

This function will calculate the mass moment of inertia of a triangular ) 
iprism given its mass, base, and height, ) 

Degin 

IBarTr i := (mTri/18) » (sqr(ta) * sqr(tb)); 

?nd; 


-unction IBarRect(var mRect,ra,rb:real):real ; 

This function will calculate the mass moment of inertia of a rectangular ) 
'prism given its mass, length, and height. ) 

<egin 

IBarRect := (mRect/12) * (sqr(ra) ♦ sqr(rb)); 
nd; 


unction WhatType: integer; 

»This function determines what shape the fluid is in. Type 1 = triangle; ») 

»type 2 = triangle and rectangle; and type 3 = one triangle and two 



•ar 

aArea, alpha, bArea :real; 

'egin 

aTriangle := 0.5 * sqr(a) » tangentTheta; 

aArea := aTriangle + (a » (b - (a * tangentTheta))); 

alpha arctan(b/a); 

if tangentTheta <) 0.0 then 

begin 

bTriangle := abs((0.5 » 5 qr(b))/tangentTheta) ; 

bArea := bTriangle ♦ abs((b » (a - (b/tangentTheta)))); 

end; 

if (angle = 0.0) then (It's a rectangle.) 

UhatType := 4 

else if (angle = 90.0) then (It's still a rectangle.) 

UhatType ;= 5 

else if theta <= alpha then 

begin 

if (aTriangle > area) then 
WhatType := 1 

else if (aArea ) area) and (bArea ) area) then 

UhatType := 2 
else UhatType := 3; 

end 
else 
begin 

if (bTriangle > area) then 

UhatType := 1 

else if (bArea > area) then 

UhatType := 6 

else UhatType := 7; 

end; 

'nd; 


Yocedure TypeOne; 

.This procedure will calculate the mass moment of inertia of a triangular ) 

iprism. ^ 

>ar 

a2, b2, h, j, cx, cyrreal; 

>egin 

h := sqrt((2 * area)/tangentTheta) » tangentTheta; 
j := sqrt((2 » area)/tangentTheta) ; 
a2 := (j/3) - (a/2); 
b2 := (h/3) - (b/2); 

CounterRotate(cx, cy, a2, b2); 

moment := IBarTri(mass,h,j) + (mass *(sqr(cx) t sqr(cy))); 

end; 


■rocedure TypeTwo; 

;This Procedure will calculate the mass moment of inertia of a volume which) 
'can be broken down into one triangular and one rectangular prism. ) 

ar 

area2, (Area of Rect) 

h, (Height of Rect) 

massTri, (Mass of triangular prism) 

massRect, (Mass of Rect. prism) 

momentTri, (Moment of inertia of triangular prism rel. to origin) 

momentRect, (Moment of inertia of Rect. prism rel. to origin) 

j, gx, gy, xbar, ybar, 

X, y, x3, y3;real; 

-egin 

h :* (area - atriangle)/a; 

j := a » tangentTheta; 

massTri := atriangle * c » rho; 

massRect := (area - atriangle) » c * rho; 

gx := a/3 - a/2; 

Qy := h + j/3 - b/2; 



Coun terRotate(xbar ,ybar ,gx .gy) ; 

momentTri := IBarTri<ma5sTri ,a,j) + (massTri » (sqr(xbar) ♦ sqr(ybar))); 

y n * 

y •= h/2 - b/2; 

CounterRotate<x3,y3,x ,y) ; 

momentRect := IBarRect(ma55Rect,a,h) ♦ (massRect * (sqr(x3) ♦ sqr(y3))); 
moment := momentTri + momentRect; 

^nd; 


Procedure TypeThree; 

iThis Procedure will calculate the mass moment o-f inertia of a volume which) 
(must be broken down into one triangular and two rectangular prisms. ) 

var 

hj Jj kj ffii 

massTri, fmass of triangle) 

massRectl, (mass of rectangle below triangle) 

massRect2, (mass of rectangle next to triangle) 

momentTri , 

momentRecU , 

momentRect2, 

tri, (the triangular area not filled with fluid) 

xl, x2, x3, x4, x5, x6, 
yl, y2, y3, y4, y5, y6, 
mcmentl, moment2:real ; 
begin 

tri ;= <a * b) - area; 
j := sqrt((2 * tri)/tangentTheta); 

h := j « tangentTheta; 

k := a - j; 

m := b - h; 

massTri ;=0.5*h*j*c* rho; 
massRectl ;= a * m * c * rho; 

massRect2 := h * k » c * rho; 
x2 != k ♦ j/3 - a/2; 
y2 := b/2 - (2 * h/3); 

CounterRotatefxl ,yl ,x2,y2) ; 

motnentTri ;= IBarTrifmassTri ,h,j) + (massTri » (sqr(xl) + sqr(yl))); 

x5 := 0; 

y5 := m/2 - h/2; 

CounterRotate(xd,y6,x5,y5) ; 

momentRectl := IBarRect(massRectl ,a,m)+(massRectl * <sqr(x6) ♦ sqr(y6>)); 

x4 := k/2 - a/2; 
y4 ;= b/2 - h/2; 

CounterRotate(x3,y3,x4,y4) ; 

momentRect2 ;= IBarRect(roassRect2,h,k)+(massRect2 « (sqr(x3) ♦ sqr(y3))); 

moment := momentTri ♦ momentRectl + momentRect2; 

end; 

Procedure TypeFour; 
var 

c4:real ; 
begin 

c4 ;= percent » b; 

moment ;= IBarRect(mass,a,c4) ♦ (mass » sqr((b/2) * (percent - 1)>); 

end; 


Procedure TypeFive; 
var 

c5:real ; 
begin 

c5 := percent » a; 

moment := IBarRect(mass,c5,b) + (mass * (sqr((a/2) * (percent - 1)))); 

end; 



Procedure TypeSix; 

{This Procedure will calculate the mass moment ot inertia of a volume which) 
(can be broken down into one triangular and one rectangular prism. ) 

uar 

area2, (Area of Rect) 

areaTri angle, 

h, (Height of Rect) 

massTri, {Mass of triangular prism) 

massRect, (Mass of Rect. prism) 

mcmentTri, (Moment of inertia of triangular prism rel. to origin) 

momentRect, (Moment of inertia of Rect. prism rel. to origin) 

j, gx, gy, xbar, ybar, 

X, y, x3, y3;real ; 
begin 

areaTriangle := 0.5 « sqr(b) * tangentBeta; 
h ;= (area - areaTriangle)/b; 

J b • tangentBeta; 

massTri := areaTriangle » c » rho; 

massRect := (area - areaTriangle) » c » rho; 

gx := h ♦ j/3 - a/2; 

gy := h/3 - b/2; 

CounterRotate(xbar ,ybar ,gx ,gy) ; 

momentTri ;= IBarTri(massTri ,b,j) + (massTri » (sqr(xbar) ♦ sqr(ybar))); 

X := (h/2) - (a/2); 
y := 0.0; 

CounterRotate(x3,y3,x ,y) ; 

momentRect := iBarRect(massRect,b,h) + (massRect » (sqr(x3) ♦ sqr(y3))); 
moment := momentTri ♦ momentRect; 

end; 


Procedure TypeSeven; 

(This Procedure will calculate the mass moment of inertia of a volume which) 

(must be broken down into one triangular and two rectangular prisms. ) 

var 

h, J, k, m, 

massTri, (mass of triangle) 

massRecU, (mass of rectangle below triangle) 

massRect2, (mass of rectangle next to triangle) 

momentTri , 

momentRect! , 

momentRect2, 

tri, (the triangular area not filled with fluid) 

xl, x2, x3, x4, x5, x6, 

yl, y2, y3, y4, y5, y6, 

moment!, moment2:rea1 ; 
begin 

tri := (a * b) - area; 
j := sqrt((2 » tri)/tangentBeta) ; 

h := j * tangentBeta; 

k := b - j; 

m := a - h; 

massTri := 0.5 * h » j » c * rho; 

massRect! := b » m « c > rho; 

massRect2 := h » k ♦ c » rho; 

x2 := a/2 - (2 * h/3); 

y2 := h/2 - (2 * j/3); 

CounterRotate(x! ,yl ,x2,y2) ; 

momentTri ;= IBarTri(massTri ,h,j)+(massTri * (sqr(x!) ♦ sqr(y!))); 
x5 := m/2 - a/2; 
y5 := 0; 

CounterRotate(x^,y6,x5,y5) ; 

momentRect! := IBarRect(massRectl ,b,m)+(massRect! » (sqr(x6) ♦ sqr(y^))); 

x4 := a/2 - h/2; 
y4 := k/2 - b/2; 

CounterRotate(x3,y3,x4,y4) ; 

momentRect2 := IBarRect(massRect2,h,k)+(massRect2 * (sqr(x3) + sqr(y3))); 

moment := momentTri + momentRect! + momentRect2; 



end; 


{»*»*«»*»»»»***»»*»»«**»*«»«»»»»»*»*»»»»»»«»**»**»»»»»»*»»»»»*»♦»♦»♦*»»****) 
Begin (*START hWIN PROGIWM*) 

ClearScreen; 

writeln <'Enter dimensions for rectangular tank;'); 
readln (a,b,c); 

writeln ('Enter water level as a percentage of 
readln (percent); 

while ((percent )= 1) or (percent (= 0)) do 
begin 

writeln('0.0 ( Water Level ( 1.0. Try again.'); 
readln(percent); 

end; 

writeln('Enter angle theta increment.'); 

readln( increment); 

writeln('Enter output file name.'); 

readln(outFileName); 

angle := 0.0; 

area := a * b » percent; 

mass ;= a * percent * b » c; 

rewr i te(outFi le,outPi leName) ; 

writeln(outFile,'A = ',a;8:4); 

writeln(outFile,'6 = ',b;8:4); 

writeln(outFile,'C = ',c:B;4); 

writeln(outFile, 'Percent full = ',percent:B:4>; 

writeln(outFile,'Type Theta Moment'); 

while (angle (= PO.O) do 

begin 

moment := 0.0; 

beta := ((90 - angle)/3i0)*(2 » pi); 

theta ;= (angle/360)»(2 * pi); 

tangentBeta := sin(beta)/cos(beta) ; 

tangentTheta := sin(theta)/cos( theta) ; 
form ;= UhatType; 
case (form) of 
l:TypeOne; 

2:TypeTwo; 

3:TypeThree; 

4:TypeFour; 

5:TypeFive; 

6:TypeSix; 

7;TypeSeven; 

end; 

writeln(outFile,form, ' ',angle:6:2, ' ', moment ;8;4); 

angle := angle * increment; 

end; 

close(outFile); 

End. (»END miN PROGRAM*) 



